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Abstract 

We develop an iterative method for finding solutions to the hermitian Yang-Mills 
equation on stable holomorphic vector bundles, following ideas recently developed 
by Donaldson. As illustrations, we construct numerically the hermitian Einstein 
metrics on the tangent bundle and a rank three vector bundle on P^. In addition, 
we find a hermitian Yang-Mills connection on a stable rank three vector bundle 
on the Fermat quintic. 
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1 Introduction 

The modern study of compactification of higher dimensional theories can be divided into 
two general branches. The first makes use of compactification manifolds with a good deal of 
symmetry, such as the torus, sphere, squashed spheres and so on. Such spaces have explicitly 
known metrics, allowing explicit solutions of the equations of motion, and explicit Kaluza- 
Klein reduction. These solutions have many applications, such as supergravity duals of large 
gauge theories; however their high degree of symmetry tends to be a problem in trying 
to obtain models with the level of complexity of the Standard Model or its often-postulated 
extensions. 

The second branch makes use of manifolds for which the relevant metrics are known 
to exist by general theorems, but for which explicit expressions are not known. The most 
famous examples are the Ricci-fiat Kahler metrics conjectured to exist by Calabi and proven 
to exist by Yau [1]. In 1985, it was proposed by Candelas et al [2] that compactification 
of the heterotic string on a Calabi- Yau manifold could lead to quasi-realistic theories of 
particle physics, containing grand unified extensions of the Standard Model and low energy 
supersymmetry. Since then, other metrics of this type, such as G2 holonomy metrics, have 
been used in quasi- realistic compactifications; see for example [3]. 
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Over the subsequent years, many tricks were developed to bypass the difficulties posed 
by not knowing the compactification metric. These tricks began with the algebraic geometry 
behind the theorems of Yau and Donaldson- Uhlenbeck-Yau, and gradually evolved into entire 
branches of mathematical physics, such as topological string theory and special geometry. 
To drastically oversimplify, the general picture is that certain "protected" quantities in the 
four dimensional Effective Field Theory (EFT), such as the superpotential in theories with 
four supercharges, and the prepotential in theories with eight supercharges, can be computed 
using techniques combining algebraic geometry with physical ideas. Other quantities, such 
as the Kahler potential in theories with four supercharges, cannot be computed directly. 
Since a good deal of important physics depends on the Kahler potential - precise values 
of particle masses, and the existence and stability of supersymmetry breaking vacua, this 
situation is not very satisfactory. 

Almost all present knowledge about the Kahler potential in the EFT comes from studying 
expansions around more computable limits. The best known example is the case of = 1 
compactifications which contain N = 2 subsectors, such as heterotic (2, 2) models, or type II 
on Calabi-Yau orientifolds. In these cases, there is a limit in which part of the N = 1 Kahler 
potential becomes equal to that of the related N = 2 theory, which is computable using 
special geometry. Other examples include the solvable orbifold or Gepner model limits, at 
which the entire Kahler potential is computable in principle using GET techniques. However, 
it is not clear at present how representative such results are of the general case. Even a limited 
ability to compute in the general case would allow studying this question. 

One completely general technique for addressing such problems is to compute the Ricci- 
flat metrics and related quantities numerically. Numerical methods are unavoidable in other 
areas of physics, beginning with such seemingly elementary problems as computing the spec- 
trum of the helium atom or integrating Newton's equations for the three body problem in 
celestial mechanics; it would be surprising if string theory could avoid this. To bring string 
theory closer to a possible confrontation with real data, for example from collider physics, it 
may be valuable to develop these missing parts of the theory of compactification. 

In this work, we make a start in this direction by showing two things. First, we re- 
view how to use existing mathematical techniques to numerically approximate metrics on 
Kahler manifolds, along lines recently developed by Donaldson [4]. Second, we extend these 
mathematical techniques to hermitian Yang-Mills connections. It will be clear that these 
techniques could be pushed to compute higher order terms, metrics on moduli spaces, and 
the like. A subsequent paper will explain the numerical methods in more detail and do 
some simple computations of terms in the EFT for compactification on a quintic Galabi-Yau 
3-manifold. 

Our direct inspirations are Donaldson's work [4] on numerical approximation of metrics, 
and of Wang [5] developing the corresponding mathematics for vector bundles. We can 
also cite Headrick and Wiseman [6], who made a pioneering numerical study of the K3 
metric using position-space methods. Finally, the ffist author is particularly indebted to 
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Bernie Shiffman and Steve Zelditch for teaching him the basics of asymptotic analysis on 
holomorphic hne bundles, and for advice in the early stages of this project, in particular for 
pointing out Wang's work. 

Let us briefly explain the problem and survey some of the approaches one might take 
towards it, before beginning the detailed development in section 2. Following [2], the deriva- 
tion of the matter Lagrangian in a heterotic compactification on a Calabi-Yau X carrying a 
bundle V involves the following steps: 

1. Find the Ricci-fiat metric Qij (with specified moduli) on X. 

2. Find the hermitian Yang-Mills connection Ai on V . 

3. Find the zero modes of the Dirac operator. As is standard, on a Kahler manifold 
this amounts to finding harmonic differential forms ip valued in V, i.e. solutions of 
= (9 + A)il) = {d + Ayip., where * denotes the adjoint operator. 

4. Find an orthonormal basis of forms if). 

5. Compute the integrals over X of wedge products of these forms to get the superpoten- 
tial. 

The key step for us is (4). Existing methods for computing the superpotential, such as 
[7, 8, 9], accomplish step (5) without needing the results of (1) and (2), by using unnormalized 
zero modes. This leads to a superpotential defined in terms of fields whose kinetic term is 
obtained from "some" unknown Kahler potential. To do better, we must either derive 
normalized zero modes in (4) for use in (5), or else take the zero modes used in (5) and 
compute their normalizations using the explicit metric from (1). 

There seems to be no way of doing this without some knowledge of the Ricci-fiat metric 
and thus the first step is to choose some approximation scheme for this metric. One's first 
thought might be to follow standard practice in numerical relativity, as done in [6], and 
introduce a six dimensional lattice which is a discrete approximation to the manifold X; in 
other words a position space approach. Taking the Kahler potential K as the basic dynamical 
variable, Einstein's equations reduce to the complex Monge-Ampere equation 

deX^ddK) =nAn (1) 

which can be solved by relaxation methods. One would then need to find similar lattice 
approximations for the connection on V and the zero modes. 

An alternative approach, introduced by Donaldson [4], is to use the natural embedding 
of X into P^^i provided by the sections of an ample line bundle (we will explain this 
in detail below). We then take as a candidate approximating metric on X the pull-back of a 
Fubini-Study metric on P^~^. Such a metric is defined by an A^ x A^ hermitian matrix. By 
suitably choosing this matrix we can try to make the associated Fubini-Study metric restrict 
to X in such a way that it gives a good approximation to the Ricci-fiat metric on X. 
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A major advantage of this approach is that it avoids the comphcations and arbitrari- 
ness involved in choosing an exphcit discretization of X; rather the entire approximation 
scheme follows from a single parameter k, the scale of the first Chern class of L. Subsequent 
mathematical development reveals more structure which can be used to our advantage. For 
example, a very natural approximation to the Ricci-flat metric, which becomes exact as 
A; ^ oo, is the so-called "balanced" metric. In a sense, to be described below, this is the 
metric for which the embedding of X into P^~i has its center of mass at the "origin". It 
also satisfies a simple fixed point condition which can be used for relaxation, solving step 

(!)• 

Another advantage, which is key for the present application, is that Donaldson's method 
can be naturally extended to study holomorphic vector bundles on X. There is a standard 
relation between holomorphic connections and hermitian metrics, which we review in sec- 
tion 2, in which step (2) of the above prescription is turned into the problem of finding 
a hermitian-Einstein metric on a vector bundle. For illustrative purposes we will explicitly 
study hermitian-Einstein metrics on two spaces: complex projective space P" and the Fermat 
quintic threefold. 

The organization of the paper is as follows. In Section 2 we provide an overview of the 
geometric background needed for our construction, in particular we will describe Donaldson's 
approach for getting metrics of constant scalar curvature. In Section 3 we explain a numerical 
approximation to the hermitian Einstein metric on a holomorphic vector bundle by a simple 
adaptation of Donaldson's scheme, building on mathematical work of Wang. In section 4 
we focus on several explicit examples. Here we describe some of our numerical methods and 
results in detail. By design we are also able to test our approximation scheme for TF'^{k), 
where TP^ is the holomorphic tangent bundle of P^, since in this case one has an analytic 
solution. 



2 Metrics of constant scalar curvature 

We follow the plan outlined in the introduction, beginning with step (1). Let X be an n- 
dimensional complex Kahler "compactification manifold." Since we are not assuming it is 
a valid string theory background, we can generalize the discussion to arbitrary n and first 
Chern class Ci{X). 

The basic example we have in mind is the complex projective space P", parameterized 
by the standard homogenous coordinates {Zj}^^Q, up to the identification {Zi} = {XZi} for 
A G C*. The Kahler potential 

KFs = k ^og(J2\Z^\'^ (2) 

defines the Fubini-Study metric on P", with SU{n + 1) symmetry Ql^j-, and g G 

U{n + 1). The parameter k controls the Kahler class uj = ddKps- 
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Our other general example is the hypersurface defined by the vanishing of a degree d 
polynomial in P": 

/(Z)= ^c--''Z,,...Z,,. (3) 

h---id 

For n = A and c? = 5 we get a quintic threefold Q. Its complex structure is determined 
by the 126 parameters c*^" *^, modulo the action of GL{5,C) on the Zj's, which leaves 101 
parameters. The generic member of this family is smooth, and has 6^'^ = 1. Therefore the 
Kahler class is determined by a single real number. A simple one parameter family of Kahler 
metrics on Q is obtained by pulling back the Fubini-Study metric on P^, or equivalently 
interpreting Eq. (2) as a Kahler potential on Q. Of course this will not be Ricci-flat. 



2.1 Approximating Ricci-flat metrics by projective embedding 

We now want to find a larger space of Kahler metrics in which to find a better approximation 
to the Ricci fiat metric. One simple generalization of Eq. (2) can be obtained by choosing 
an {n + 1) x (n + 1) hermitian matrix h^^ , and writing 



K^, = k\og\ Y^h^^Za, 

J,j=o 



(4) 



Of course, by making a linear redefinition of coordinates, we could turn this back into Eq. 
(2), but doing so would modify Eq. (3). Rather, by fixing Eq. (3), this way we get an 
(n + 1) ^-parameter family of Kahler potentials. 

Another way to think about this definition is to make the linear redefinition taking h to 
the identity. In this case, the parameters we are varying to control the metric are the extra 
25 parameters in Eq. (3) determining a specific embedding of Q into P^. While all of the 
embeddings are equivalent under a GL{5,C) action, once we use the metric, we break this 
to U{5)] thus the set of metrics we can obtain this way is parameterized by a GL[5, C)/U{5) 
homogeneous space. 

A simple generalization to get more parameters could be motivated by noticing that Eq. 
(2) is also equal to 



Kfsj 



log E 



(5) 



and generalizing this to 



log j: 



,ik=0 



(6) 



h,k 



log 



E 



7- 



,ik,jl,--- ,jk=0 



]k 



(7) 



which can again be interpreted as a Kahler potential on Q. In simple terms, we are using 
higher degree polynomials as basis functions. Now we have an (n + l)^'^-parameter family 
of metrics, and by taking k large we can imagine finding an arbitrarily good approximating 
metric within this class. 

One way to find the best approximation to the Ricci-fiat metric on Q would be to write 
Eq. (1) directly in these variables. Note that the holomorphic (n, 0)-form Vt is known 
explicitly. For example, in the coordinate patch where 7^ we can choose the local 
coordinates Wi = Zi/Z^, in terms of which 

^ dwi dw2 dw3 
df/dwi 

and thus one can write the volume form for the Ricci fiat metric explicitly, 

dvolx = nACt (8) 

without solving any equations. One might then substitute Eq. (7) into Eq. (1), evaluate 
this at a set of points Pi, and solve the resulting system of nonlinear equations. These are 
rather complicated, however, and furthermore we have introduced arbitrariness in the choice 
of the Pi. Now this arbitrariness can have its uses, for example we might use it to place 
more points in regions of large curvature. On the other hand, it means that the results will 
not have simple mathematical or physical interpretations, except in the limit in which the 
number of points is so large that we can ignore the discretization.^ Before investing a lot of 
effort into their study, we should try to improve on this point. 

2.2 Balanced metrics 

There is a pretty construction that goes back to [10] which provides a more natural approx- 
imating metric, and a numerical scheme which is guaranteed to converge to it. 

First, we can systematize the construction which led to Eq. (7), by noting that the basis 
functions are products of degree k holomorphic times degree k antiholomorphic monomials. 
Let the number of independent holomorphic degree k monomials be iV+l; this is the binomial 
coefficient {"'^^) for P", and we will give it for Q later. 

Let us phrase this construction in a way which can be used for an arbitrary manifold X. 
We choose a holomorphic line bundle C over X, with global sections. Denote a complete 
basis of these as Sa, where 1 < a < N, and consider the map 

ik-.X^ pA^-l ,,(Zo, . . . , Z„) = S2iZ), . . . , SNiZ)). 

The geometric picture is that each point in our original manifold X (parameterized by the 
Zi) corresponds to a point in parameterized by the sections Sa- Since choosing a different 

^Or unless we can come up with a construction in which some sort of physical objects at the points pi 
enforce the equations. 
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frame for L would produce an overall rescaling — > Asq,, the overall scale is undetermined. 
Granting that S2{Z), . . . , Si\[{Z) do not vanish simultaneously, this gives us a map to 

pAf-l 

The simplest example is to embed using C = Opi{k) into F^. In this case the map is 
ikiZo, Zi) = (Zq, Zq ^Zi, Zq ^Zf, ZqZ'I ^, Z^). 

In general we want this map to be an embedding, i. e. that distinct points map to distinct 
points with non-vanishing Jacobian. In general, we can appeal to the Kodaira embedding 
theorem, which asserts that for positive C this will be true for all for some k > k^. For 
non-singular quintics, this is true for OM{k) for all A; > 1. As a point of language, the pair 
of a manifold X with a positive line bundle C is referred to as a polarized manifold {X, C); 
the condition that this construction provides an embedding for some k is that C is ample. 

Now, we consider our family Eq. (7) of candidate Kahler potentials, and rewrite them as 

J^, = log l^r^s^s^ 
or simply 

Kh = \og\\s\\l (9) 

for short, where Sa plays the role of a degree k monomial. We now have an A^^-parameter 
family of Kahler potentials, and will seek a good approximating metric in this family. Just 
as before, this amounts to using the pull-back of a Fubini-Study metric from as our 

trial metric. 

Mathematically, the simplest interpretation of Eq. (9) is that it defines a hermitian 
metric on the line bundle C = OM{k)- This is a sesquilinear map from C ® C io smooth 
functions C°°(X), here defined by 

(s, s') = e-^'^ -s-s' - ^ ' ^ 



The point is that a change of frame, which acts on our explicit sections as Sa — > Asq,, cancels 
out of this expression.^ 

This metric allow us to define an inner product between the global sections: 

Hap = (splsa) = i / ^f^dvolx. (10) 
Jx IpIU 

This is the "physical" inner product in a sense we will explain further below. Note that it 
depends on /i in a nonlinear way, since h appears in the denominator. 

possibly more familiar physics use of this is in = 1 supcrgravity: taking K —K and s W, one 
gets the standard expression for the gravitino mass e^\W\'^. In an example such as the flux superpotential, 
in which is a sum of various terms Sa with constant coefficients, Eq. (9) also applies to give K. 
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Here dvolx is a volume form on X, which has to be chosen. If X is Calabi-Yau, it is 
simplest to use Eq. (8) to define dvolx- If X is not Calabi-Yau, the standard choice of dvolx 
is to take dvol^^ = uj^/n\, where u is the Kahler metric derived from Eq. (9). This depends 
on h as well, so the expression is even more non-linear in h. 

Thus, given h and a basis of global sections s^, we could compute the matrix of inner 
products Eq. (10). Once we have it, we could make a linear redefinition, say s = H^^^'^s, 
and go to a basis of orthonormal sections where 

= ^o.B- (11) 

On the other hand, Eq. (9) also implicitly defines a notion of orthonormal basis locally 
in the bundle, in which 

/i"^ = 5"^. (12) 

This is a priori different from Eq. (11); indeed we can freely postulate it when we write Eq. 
(9). However, if the two notions agree, 

HajS = (^ ^)a/3) 

then we can go to a basis of sections in which 

H^-p = h-~^ = (13) 

In this case, the embedding of X in P^~^ using these sections is called balanced. More 
generally, we call a polarized manifold (X, C^) balanced if such an embedding exists. 

An equivalent definition of the balanced embedding is arrived at if we consider the func- 
tion on X defined as 

p{u){x) = Y,{H~Y'M^)M^)) (14) 

or equivalently 

P{^)i.x) = ^\\Sa{.x)\\^ 
a 

where the second sum is taken over an orthonormal basis in which H = dap- X is balanced 
precisely when p{ijj){x) is the constant function. 

Many theorems have been proven about balanced manifolds. Let us first recall the 
following theorem of Donaldson (Theorem 1 in [11]): 

Theorem 2.1 Suppose the automorphism group Aut(X, £) is discrete. If {X^C^) is bal- 
anced, then the choice of basis in H^{X,C^) such that ik{C) is balanced is unique up to the 
action ofU{N) x W. 

The condition on Aut(X, £), i.e., there are no continuous symmetries, is true for the quintic 
Q. This theorem then tells us that, if a metric h exists which gives a balanced embedding, 
it is unique up to scale. 
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Given a balanced embedding, one defines the balanced metric on X as the pullback of 
the Fubini-Study metric (9): 

27r ^ , , - - 

= -^iki^Fs), (15) 

The cohomology class of the Kahler form [u}k\ = 27rci(£) G H'^{X) is independent of k. 
Using these definitions Donaldson proves that (Theorem 2 in [11]): 

Theorem 2.2 Suppose Aut(X, £) is discrete and {X^C^) is balanced for sufficiently large 
k. If the metrics Uk converge in the C°° norm to some limit Uoo as k ^ oo, then uJoo is a 
Kahler metric in the class 2tcci{C) with constant scalar curvature. 

The constant value of the scalar curvature is determined by ci{X), and in particular for 
ci{X) = the scalar curvature is zero. Thus, the balanced metrics ujk, in the large k limit, 
converge to the Ricci flat metric. 

Therefore, if we can find the unique balanced metric for a given C, it is a good candidate 
for approximating the Ricci fiat metric on X. One may ask where the complex structure and 
Kahler moduli on which this Ricci fiat metric depends, are put in. The complex structure 
enters implicitly through the basis for holomorphic sections Sa, as we will see in examples 
below. As for the Kahler class, recall that this is determined, up to scale, to be 27rci(£). 
Of course, the Ricci flatness condition is scale invariant, so the overall scale is irrelevant; 
however the point of this is that if 6^'^ > 1, then by appropriately choosing C we choose a 
particular ray in the Kahler cone. This will not be relevant for our examples here but shows 
that in principle any Ricci-flat Kahler metric could be approximated in this way. 

2.3 Finding the balanced metric 

In [12, 4] Donaldson proposes a method to determine the hermitian metric h in Eq. (9), 
which will lead to a balanced metric. He defines the "T operator" , which given a metric h 
computes the matrix H: 

Hai3 = T{h)^p = ^ I rfvolx (16) 

Now, suppose we find a fixed point of this operator, 

T{h) = h. 

Then, by a GL{N) change of basis s h''^^^s, we can bring h to the unit matrix, which 
will produce the balanced embedding. 

The simplest way to find a fixed point of an operator is to iterate it. If the operator is 
contracting, this is guaranteed to work. In our case we have the following theorem [11, 13]: 

Theorem 2.3 Suppose thatAnt{X,L) is discrete. If a balanced embedding exists then, for 
any initial Gq hermitian metric, the sequence T^{Go) converges to the balanced metric G as 
r — i> oo. 
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Thus the T operator can be used to find approximate Ricci-fiat metrics on Calabi-Yau 
manifolds, and more generally approximate constant scalar curvature Kahler metrics. In [4] 
Donaldson studies numerically explicit and K?) examples. We will discuss some additional 
examples below. 

2.4 Balanced metrics and constant scalar curvature 

In this subsection we outline the reason why the limit of a family of balanced metrics has 
constant scalar curvature. This is the content of Theorem 2.2. This will be very useful later 
on, when we generalize the T-operator to vector bundles. 

Note that the function p{uj) is independent of the choice of orthonormal basis, and remains 
unchanged if we replace /i by a constant scalar multiple. Therefore, it is an invariant of the 
Kahler form. As discussed before, the balanced condition for (X, L^) is equivalent to the 
existence of a metric ujk such that p{u}k) is a constant function on X. The asymptotic 
behavior of the "density of states" function p{ujk) as A; — > oo for fixed uj has been studied in 
[10, 14, 15, 16]. Note that for any metric 



J X 

where the coefficients can be determined using the Riemann-Roch formula. Note that ao 
is just the volume of X and 



where S{uo) is the scalar curvature of uj. We will use the following result (Prop. 6 in [11]): 
Proposition 2.4 1. piuo) has an asymptotic expansion as k oo 




(17) 




Pfc(cu) ~ Ao(cu)p + Ai(cu)r-^ + ■ ■ ■ 



where Ai{uj) are smooth functions on X defined locally by u. In particular, 



Aoiu) = 1, A,{u) = —S{u). 

ZTT 



2. The expansion holds uniformly in the C°° norm; in that for any r,N>(} 



N 



.n—i 



.n-N-1 



i=0 



for some constants Kr^N,, 
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Now assume that we are given balanced metrics uJk converging to u^. Then by the 
previous proposition 



p,K) - A;" - —S{uJk)k''-^ 
Ztx 



^ ck 



n-2 



for some constant c. Since is balanced Pk{^k) is constant: Pki'^k) = dimif'^(X, L'^)/\^, 
and we can use (17) to find that 



V Ztt 



^ ck 



n-2 



CO(X) 



or equivalently 



0{k~'] 



CO(X) 



Hence S'(ci;oo) = •S'o = constant, where Sq = y Jj^ S{uj) is the mean curvature. 



3 The hermitian Yang-Mills equations 

We are now ready to generalize the T-operator, which provided an approximation scheme 
for the constant curvature metric, to a "generalized T-operator" which can be used to find 
a solution of the Yang-Mills equations on a Calabi-Yau manifold X. 

We briefiy recall the argument that a solution of the Yang-Mills equations which preserves 
A/" = 1 supersymmetry, must be hermitian Yang-Mills. First, the supersymmetry variation 
of the gaugino has to vanish, 

where F^j, is the Yang-Mills field strength, and e is the covariantly constant spinor. 
Going to complex coordinates (i, i) and rewriting of the Clifford algebra as 

Ti dz'; Tj d^, 

this is equivalent to 

This is the particular case of the hermitian Yang- Mills equations with Tr F = 0. The 
general case replaces the last equation with 

u^Ffj = c ■ 1 

for a constant c, determined by the first Chern class. For convenience we abbreviate this 
equation below as 

^F=c.l. 
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Next we review the relation between solutions of these equations, and holomorphic bun- 
dles carrying hermitian-Einstein metrics. In physics, one defines Yang-Mills theory in terms 
of a connection on a vector bundle with a fixed metric. First, a connection on a vector bundle 
can be described in terms of a connection one-form by choosing a frame for the bundle, say 
ea{x), and defining the covariant derivative as 

Div'^ea) = {dv^)ea + v''Aleh. 

In physics, one usually takes the frame to be orthonormal, (e^, e^) = 6ab, and thus 

K^O = Ke„^;''e,) = K)*^;^ (18) 

where * is complex conjugation. 

The condition that the connection be compatible with the metric, 

d{u,v) = {Du,v) + {u,Dv), (19) 

reduces to requiring the connection one-form to be anti-hermitian. 

In mathematics, one often considers a more general frame, for which the metric is a 
hermitian matrix, 

(ea,eb) = G,b, G = G\ (21) 
Decomposing the positive definite hermitian matrix G as 

G = h^h, (22) 

we see that the math and physics conventions differ by a complex gauge transformation: 
u = hs. This complex gauge transformation leads to a different form for the connection, 
according to the standard relation 

di + A'^'^'''^\ = h{d, + A^P''y'\)h~\ (23) 

Now, equations of the form 

will be integrability conditions for the covariant derivatives. In particular, this equation has 
the general solution 

in other words the D covariant derivatives are obtained from the derivative (9 by a complex 
gauge transformation. 

Thus, we can use Eq. (23) to bring the connection to the gauge ^(™«*'*) = o, at the cost 
of losing the simple metric Eq. (18) and Eq. (20). Actually, the covariant derivative is still 
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compatible with the metric as in Eq. (19), we just have a non-trivial fiber metric G. The 
metric compatibility condition becomes 

= d{u, v) = {Bu, v) + (u, Dv) 

so 

or equivalently 

^{math) ^ G^^dG. 

Conversely, if we are given a metric G, then we can use the inverse complex gauge 
transformation to bring the connection back to the unitary form. This leads to the formula 

Afy^'^ = h{dih~^). 

Using Eq. (20), we can get the entire connection, so the metric contains the same information 
as a connection satisfying F*^'''^-' = F*^^''') = 0. Thus we can rephrase the final equation on 
F'^^'^\ as a condition on the metric. It is simplest to write this in the "mathematical" gauge 
^{math) = 0, in which it is 

c ■ 1 = uj'^^Ffj = J'd-jA^'^'''^\ = J'd-j {G~^diG) . (24) 

A metric G satisfying this equation is a "hermitian-Einstein" metric. It is simply related to 
a hermitian Yang-Mills connection as above. 

Finally, using the complex gauge transformation above, the standard physical inner prod- 
uct 

{u\v) = [ {u*)av'' (25) 
Jx 

is equal to the natural inner product generalizing Eq. (10), 

{u\v) = {h s\h t) 

= z I G^isU^dYolx 
Jx 

3.1 Embedding vector bundles 

We now want to represent the hermitian metric in the same way as we did for line 
bundles, by introducing a complete basis of sections. Now an irreducible bundle E with 
Ci = 0, and thus of interest for string compactification, will not have global sections. What 
we do instead is to make the same construction for E{k) = E ® C^, which will have global 
sections. We can again think of these sections as a basis of polynomials approximating 
functions on which to base our numerical scheme. 

Thus, consider a rank r vector bundle E, and suppose that E(k) has global sections. 
Choosing a local frame as above, a basis for these will be an by r matrix z^. This is 
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(26) 
(27) 



defined up to a GL{N) change of basis, and up to a GL{r) change of frame. After making 
these identifications, such a matrix z defines a point in the Grassmannian G{r, N) of r planes 
in C^. 

Given a metric G^j, on the fibers of E{k), we can define the matrix of inner products 

as above. Such a metric could be obtained by multiplying a metric G^^^ on E, by one on 
as defined earlier. Or, it might simply be an r x r hermitian matrix of functions (in each 
local frame) with appropriate transformation properties. 

Now there is a natural set of metrics on E{k) generalizing Eq. (9), again parameterized 
by an X matrix, defined by 

^G-Y=g-^z^^{z%, 

where the dagger is hermitian conjugation. Again, the approach will be to find a natural 
metric in this class which is a good approximation to the hermitian-Einstein metric. This 
will lead to a hermitian Yang-Mills connection on E{k). But this is simply related to the 
hermitian Yang-Mills connection on E, because twisting by only modifies the trace part 
of the field strength. 

3.2 Generalized T-operator 

We will now turn to a proposal for a generalized T-operator, which produces the hermitian- 
Einstein metric on a stable vector bundle. To begin with we use results by Wang about 
balanced metrics on such bundles [5]. 

We consider again a polarized n dimensional manifold (X, C) and an irreducible holo- 
morphic vector bundle E of rank r on X. Then by Kodaira embedding we know that for k 
sufficiently large, a basis z'^ of the global sections of E{k) will give rise to an embedding 

X^^G(r,A^). 

Now Wang proves the following: 

Theorem 3.1 E is Gieseker stable iff' there is an integer ko such that for k > ko, the kth 
embedding given as above can be moved to a balanced place, i.e., there is a g E SL{N,C) 
which is unique up to left translation by SU (N) such that: 

If r 

— / z{zhy^z^ dV = —InxN- 

^ Jg X 

We call the equation above the "balance equation." In the case that E is & line bundle, 
this definition reduces to that of a balanced embedding in P^~^. 
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Now, let h he a hermitian metric on C and if be a hermitian metric on E, and fix the 

Kahler form on X to be = — Ric{h). Let vol denote tlie volume of {X,uj). Suppose 

Si, . . . , Sn is an orthonormal basis of H^{X, E{k)) with respect to the induced -metric 
(. , .). The Szego kernel Bk is a generalization of the function p{uj) defined in Eq. (14). It is 
defined as the fiberwise homomorphism 

N 

Bk{x) = S,{x))S,{x) : ^ E^. 

1=1 

This expression is independent of the choice of orthonormal basis. 

Now the local form of Theorem 3.1 can be stated as follows (Corollary 1.1 of [5]): 

Theorem 3.2 E is Gieseker stable iff there is an integer ko such that for any k > k^, we 
can find a metric H^^\ which we will call the balanced metric on E(k), such that the Szego 
kernel satisfies the equation 

where Ie is the identity bundle morphism and x(^) is the Hilbert polynomial of E with respect 
to the polarization L. 

The theorem tells us that if E is Gieseker stable then for large k there is a balanced 
metric B.'^^^ on E{k). Hence we will have a sequence of hermitian metrics Hf^ := H^''^ ® h~'' 
on E. The importance of the balanced metric H^'^^ for physical applications follows from the 
following theorem: 

Theorem 3.3 Suppose E is Gieseker stable. If — > H^o in the C°° norm as k ^ oo, then 
the metric ifoo solves the "weak hermitian- Einstein equation", 

±M,.„.,4si.,I.^{^.l).I. (28) 

where /\F(^E.Hao) ^■^ the contraction of the curvature form of E with uj, S{uj) is the scalar 
curvature of X and s := y S^ . Conversely, suppose there is a hermitian metric H^o 
solving this equation, then — > in norm for any r. 

To prove (28) one can work along the same lines as in the proof of Theorem 2.2, using 
Catlin's and Wang's results for the expansion of the Szego kernel. 

Proposition 3.4 1. For fixed hermitian metrics H and h on E and Ox(l) respectively, 
there is an asymptotic expansion as k oo 

BkiH, h) ~ Aq{H, h)k'' + Ai{H, h)k''~^ + ■ ■ ■ , 

where Ai{H, h) eT{ End E) are smooth sections defined locally by H . In particular, 

AoiH, h) = Ie. Ai{H, h) = ^ A F{E, Ric{h)) + ]-S{u) ■ Ie 
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2. The expansion holds uniformly in the C°° norm; in the sense that for any r, N > 

N 

for some constants Kr^N,H,h- 

Now we can repeat the steps of the argument outhned in Section 2.4. Under the assump- 
tion that Hk Hao in we find that for r > 

\\Bk{H,) - lEk^ - ^ f\F{E,mc{h)) + \s{uj) ■ lEk'^-^Wcr ^ Ck^~' 

for some fixed constant C. By assumption H^'^'^ is balanced, hence Bk{Hk) = xik)/rVlE- 
This imphes that 

3.2.1 Generalized T-operator 

Using the strong analogy between the construction of metrics with constant Kahler curvature 
and metrics on stable bundles which obey the hermitian-Einstein equation, we propose the 
following generalized T-operator: 

T{G) = — / z{z^G-'z)-^z^ dV, (29) 
Vr Jx 

where as before, 2; is an by r matrix of holomorphic sections of E. 
The relevance of this proposal follows from the following conjecture: 

Conjecture 3.5 // a balanced embedding i: X ^ G{r,N) exists, then for every hermitian 
N X N matrix G, the sequence T^'{G) converges to a fixed point Gq as r ^ 00. Using an 
orthonormal basis with respect to Gq, the embedding is balanced, and as outlined above, it 
provides an approximate solution for the corresponding hermitian-Einstein equation. 

This conjecture may require additional technical assumptions, such as the earlier one of 
Aut(X, E) being discrete. We have not attempted to prove it, but would hope that this can 
be done along the same lines as [11, 13]. 

In the following section we will numerically test the conjecture for several stable vector 
bundles on P^, and on the Fermat quintic in P^, and find that it works for these cases. 
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4 Examples 



4.1 Hermite-Einstein metric on the tangent bundle of 

Let P" be the complex projective space of dimension n, and {^j} -=0 homogeneous coordi- 
nates. We will work on the open set Zq ^ and chose the local inhomogeneous coordinates 
Wi = Zi/Zo. The Fubini-Study metric on P" 

is the unique maximally symmetric metric, with its group of Killing symmetries isomorphic 
to U{n + 1). In addition, this metric is Einstein, that is its Ricci tensor is proportional to 
the metric itself. Therefore its associated curvature tensor obeys the hermitian Yang-Mills 
equation. The Donaldson-Uhlenbeck-Yau theorem then implies that the tangent bundle of 
P", TP", is a rank n stable bundle on P".^ It follow from this that the balanced metric on 
the bundle TP"- must be the Fubini-Study metric. 

To describe the tangent bundle TP" we use the Euler sequence 

— , Op,, — , C»pn(l)®("+i) — > TP" — > 0. (30) 

Here Opn(l) denotes the hyperplane line bundle. After twisting the sequence by Cpn(/c) and 
taking the cohomology we find the short exact sequence (SES) 

— > i7°(P",Opn(A;)) — > i7°(P",Opn(A; + 1)®("+^)) — > iJ°(P", TP"(A;)) — > 0. (31) 

This gives an explicit description for iJ''(P", TP"(A;)), which for sufficiently large k gives the 
embedding 

P" ^ G(n, W) (32) 

where W = if°(P", TP"(/c))*, and G(n, W) is the Grassmannian of n-planes in W. 

Based on (31), we choose to describe the global sections of TP"(/c) by an n + 1 vector 

(Mo, . . . , MO 

where {Mj}"^^ are arbitrary monomials of degree -|- 1 in the homogeneous coordinates Zi, 
while Mq is any degree k + 1 monomial which does not contain an Zq. 

Now we show how to construct the embedding (31) for any A; > 0. We start by choosing 
a frame {ej}"^o vector bundle 0{k + 1)®("+^). This amounts to choosing a section for 

every one the n + 1 summands. For simplicity we chose the same section in every summand. 
The Euler sequence (30) imposes the condition 

n 

^ ZiCi = 0. 
■^The stability of TP" also has purely algebraic proof. 

17 



Locally this gives a frame for TP", if we solve for 

i=i ° i=i 

Expanding the global sections of TP"(A;) in the local frame {cj}"^^ gives an n x dim(iy) 
matrix, which is the explicit realization of our embedding [17]. 

To illustrate the procedure consider TP^(O). Cp2(l) has 3 global sections: Zq,Zi,Z2- 
Choosing Zq to be the local frame in every summand of (9p2(l)®^, and discarding the global 
section Zq from the first (9p2(l), we find the matrix 

/ —vol —W1W2 1 uii 
\ —W1W2 —wl 1 Wi W2 

For an initial hermitian metric Go on the vector space W = /7°(P", TP"(/c))*, our gener- 
alized T-operator (29) gives the iterations 

We tested the converges of the T-map starting with Go = / in the case n = 2 for 
/c = 1, . . . ,5. In all cases we converged to a given Goo for less than 10 iterations, with a 
precision of 0.1%. 

The balanced metric H^'^^ on the vector bundle TP"(/c) induced by Goo is given by 

H^'^^ = {z^G^h)-\ (34) 

Let h be Fubini-Study metric on the hyperplane bundle Opn(l), that is the metric with 
constant scalar curvature. Then the metric 

is the balanced metric on TP". Our numerical computations show that this is indeed the 
Fubini-Study metric on TP", as explained earlier. The numerical agreement is within 0.5%. 
This provides the first non-trivial test of our conjecture. 

4.2 A stable rank 3 bundle over 

In this section we test our generalized T-operator on a rank 3 vector bundle V* over P^. We 
first consider its dual V, defined by four linearly independent global sections {nii} of (9p2(2) 
through the SES 

-0®^^Op2(2) -0. 

This bundle has moduli, which are implicitly determined by the choice of the sections {rrii}. 
Before choosing these, let us check stability, which does not depend on the specifics of this 
choice. 
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To check stability, we have to ensure that neither V nor A^V have destabihzing hne 
bundles. Using the canonical isomorphism 

A^y = det 1/ ® V* 

we find the slopes 

I2{V) = -2/3, ;u(aV) = -4/3. 

Since PiciF"^) = Z, all line bundles are of the form Op2{p) for some p. Hence it is sufficient 
to show that 

H''{F\V)=0, if°(p2,AV(l)) = 0. 

The first fact follows from the defining sequence of V, if we assume that {rrii} are linearly 
independent. To prove the second statement we use 

H^{¥^, aV(1)) = iJ°(P^ V*{-1)) = if2(P^ V{-2)y. 

Again, this statement follows easily from the defining sequence of V. Finally, stability of V 
implies stability for V*. 

We will now compute the hermitian Yang-Mills connection on V* using our generalized 
T-operator. First observe that V*{k) fits into the short exact sequence 

^ Or2{k - 2) ^ Op2{kf'^ V*{k) 0. (35) 

This leads to another SES 

^ H\F\ 0^2{k - 2)) ^ 0^2{k)®^) ^ /f0(p2, V\k)) ^ 0. 

We can use this expression for an explicit parameterization of if'^(P^, V*{k)). 

For concreteness let us choose to be four global sections {mj}^^^ defining V to be 

Z1Z2, ZqZi, ZqZ2, Zq. 

Now we choose a frame {cj} for Or>2[k)®^. The defining equation (35) of V*{k) imposes the 
condition Y^^rriiii = 0, and gives a frame for V*. Locally we can solve for Cq, and working 
in inhomogeneous coordinates Wi we find that 

1 1 1 

Co = ei 62 63. 

W2 Wi W1W2 

Expanding the global sections of if°(P^, V*{k)) in the frame {e^jf^^ gives a matrix, which is 
the embedding map. 

We studied the convergence of our generalized T-operator numerically for k = 2,3 and 4, 
and found that convergence was achieved for less than 10 iterations. As before, the metric 
on V*{k) is given by 

hO'^ = {z^G-Jz)-\ (36) 
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while the corresponding metric on V* is 



where h is again the Fubini-Study metric on Cp2(l). 

Since in this case the balanced metric on V* is not a priori known, one needs a different 
approach, than used in the previous section for TP^, to test how close is the approximate 
balanced metric to satisfying the hermitian Yang-Mills equation. But this quite easy to do 
numerically once the balanced metric Goo is known, as all we need to do is to check how 
close we are to satisfying Eq. (24). In all cases considered Eq. (24) was satisfied to within 
1% accuracy. 

4.3 A rank 3 bundle on the Fermat quintic 

In this section we turn to a much more complicated case than our previous examples, that 
of a stable rank 3 bundle on the Fermat quintic Q in P^: 



Testing our generalized T-operator in this case necessitates knowledge of the Ricci fiat 
metric on the Fermat quintic. For this we use Donaldson's original T-operator Eq. (16), 
which we turn to first. 

4.3.1 Ricci fiat metric on Fermat quintic 

We consider the embedding of Q given by the complete linear system of cubics, H^{Q, Cq(3)), 
whose complex projectivization is isomorphic to P^''. The balanced metric will be the re- 
striction of a Fubini-Study metric on P^^. An indirect test that this has indeed vanishing 
Ricci curvature is included in the next section. 

In order to do practical calculations with Donaldson's T-map, we have to perform the 
integrals on Q numerically. We introduce a discrete approximation to the Calabi-Yau volume 
form d^Q = A fi, defining it by a weighted set of M points {xa\^=i G Q, with masses Va'- 



This numerical measure gives an accurate approximation to the analytical one for large M. In 
our computations we build 10 different samples of 100,000 points, which we use independently 
to iterate the T-map until convergence is reached, i.e., the sequence {r''(Go)}r=o obeys 



Q: 



zl + zl + --- + zl = 0. 



(37) 




(38) 



T'+\Go)-T'{Go)\\<e. 



In our simulations the fixed point of this discrete version of the T-map was reached after 
15-20 iterations. Each weighted point set gave rise to a convergent sequence. The 10 different 
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hermitian forms {G^}g2=i approximating the balanced metric in FH'^^Q, Oq{3)) agree up to 



max 



0.9%, 



(39) 



where is the average matrix of the ten different outputs and o"(G'^) is the standard 

deviation matrix. The ratio (G^) | is computed entry by entry, and the maximum is 

taken over all entries. We used the average (G^) as approximation for the balanced metric 
on H^{Q,Oq{3)). 





Figure 1: The shape of the rational curve for the balanced and non-balanced metrics. 

To get a visual picture of the geometry implicit in the construction, we consider the 
rational curve t: Q, defined locally by the parametrization 

( 1, -1, t, 0, -t ) (40) 

with t G C U oo. Take the sections Zf + Z| and Zq from H^{Q, Oq{3)), and consider the 
function s = {Zf + Z|)/Zq = wf + wf. In Fig. 1 we consider the real function \s\q restricted 
to the rational curve (40), where we take the stereographic projection of the complex t-plane 
and for a given t we plot |s(t)||. in the radial direction. For the balanced metric (G^) the 
deviation from being spherical is small. For comparison we also show the same plot for the 
case of a generic non-balanced hermitian form G with random entries. 
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4.3.2 Solution of the hermitian Yang-Mills equation 

In this section we use the generalized T-operator to produce a hermitian Yang-Mills connec- 
tion on a rank three stable vector bundle V on the Fermat quintic Q. We also implicitly test 
that the previously obtained balanced metric on Q indeed has vanishing Ricci curvature. 
We define the rank three bundle V by the following SES 

— -0. 

/3 is given by four generic global sections of Oq{1), which do not intersect on Q, hence V is 
indeed a vector bundle. In addition, the first Chern class of V is ci{V) = H, hence V is not 
a simple twist of the tangent bundle of Q. That fact that V is stable was proved in [18]. 
Once again, we use 

^ Oq(A; - 1) — 0^\k) V{k) — 0, 

and its associated long exact sequence in cohomology 

H%Q, Ogik - 1)) H\Q, 0^\k)) ^ H\Q, V{k)) 0, 



to derive a frame for V and an explicit parameterization for the global sections. We choose 



(3 = {Zq, . . . , Zs). Using the frame {ejjf^Q for Oq^, we also get a frame for V with the 



relation 

3 

Co = - y^^Wjij. 

1=1 

In this paper we restrict to the case = 1 for which dimH^{Q,V{l)) = 19. The 
coordinate matrix 



1 . . . W4 










-WIW2 


-W1W3 


—W1W4 





1 . . . W4 





-WIW2 


-wl 


-W3W2 


—WiW2 








1 . . . 


-WIW3 


-W2W3 


-wl 


— WaWz 



z{w) =1 1 . . . W4 —W1W2 —W2 —wsW2 —W4W2 I (41) 
\ 1 . . .W4 —W1W3 —W2 

gives the embedding into the Grassmannian Q ^ ^(3, 19). 

Using the integration techniques described in the previous section, we iterate the general- 
ized T-operator. We reach the fixed point after 12-15 iterations for several different samples 
of weighted points which approximate the analytical measure, allowing us to estimate the 
balanced metric for H^{Q, ^(1)) with an error of 1.1%.^ 

The metric on V{1) is given by 

H = iz^G^h)-\ (42) 



"'We estimate the errors using (39). 
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To test the accuracy of this metric we evaluate the right hand side of the hermitian Yang- 
Mills equations (24). We find the mean to be 

i^^F^) = ^ (^^^>.j)rf/in ^ 1.31 X I3x3 

with I3x3 the 3x3 identity matrix. In our conventions the theoretical value of the constant 
is 4/3. The standard deviation of the individual matrix elements is 



cr(co'*-'-Fjj) = max 

where the square-root and the square are performed entry by entry. Therefore, oj^^Ffj is a 
global constant on Q times the identity, within an error of 0.15/1.31 ~ 11%. This implies 
that the hermitian Yang-Mills equation (24) is satisfied with this accuracy. 

Testing the hermitian Yang-Mills equation provides an implicit test of Ricci flatness, since 
it is precisely the Ricci fiat metric that is needed in the hermitian Yang-Mills equation. If we 
had gotten this metric wrong, then we would have had no chance of satisfying the hermitian 
Yang- Mills equation. 




Figure 2: The probability density in the radial direction on the rational curve. 

Finally, to visualize the construction, in Fig. 2 we took the rational curve defined in (40), 
and we plotted the function \^\%, where 

^> = {wiWijei + {w2Wi)e2 + (^3^4)63 

and G is the balanced metric we obtained. If we interpret \1/ as a wave-function, then Fig. 2 
exhibits the probability density (^'I^E') in the radial direction, restricted to the rational curve 
(40). 
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